Novel Geometrical Models of Relativistic Stars. 
I. The General Scheme 
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In a series of articles we describe a novel class of geometrical models of relativistic stars. Our 
approach to the static spherically symmetric solutions of Einstein equations is based on a careful 
physical analysis of radial gauge conditions. It brings us to a two parameter family of relativistic 
stars without stiff functional dependence between the stelar radius and stelar mass. It turns out 
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I. INTRODUCTION 

Today the theory of relativistic stars is a well developed 
branch of relativistic physics with many observational 
confirmations. It is based on general relativity (GR), 
as a relativistic theory of gravity, on quantum statisti- 
cal physics, as a theory of many-particle systems, and on 
the latest achievements of the standard model (SM), as 
a modern theory of matter constituents. (See the books 
[1] and the large amount of references therein.) 

The role of quantum statistics of the Fermi gas in the- 
ory of neutron stars has been clarified in the pioneering 
articles by Chandrasekhar and Landau [2] on the non- 
rclativistic ground of Newton gravity. For this purpose 
was used the analogy with the theory of white dwarfs [3] . 

The beginning of the relativistic stelar theory can be 
found in the pioneering articles by Schwarzschild [4] , Tol- 
men, and Openheimer&Volkov (TOV) [5], an important 
further developments - in [6] . 

After the appearance of these articles the relativistic 
theory of gravity of spherically symmetric static stars is 
widely accepted as a well established issue. The further 
developments are related with various considerations of 
the physics of stelar matter and with a search of more re- 
alistic equations of state (EOS) of this matter for different 
types of stars. This line of investigation is continuously 
followed up to now, see [7] and the references therein. 

In spite of general success of the relativistic theory of 
stars, at present it can not be considered as a complete 
established scientific area in its final form. (See, for ex- 
ample, the recent review articles [8] and the references 
therein.) 

Some difficulties in the explanation of the properties 
of very dense stars, like neutron stars or eventual quark 
ones, are suspected [9] and need a proper explanation. 
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The currently used approach of continuous modification 
of EOS already brought us to EOS, which are extension 
of known physical laws onto domain where we have no ex- 
perimental information. Therefore some additional, more 
or less arbitrary assumptions are needed. 

To some extend, a similar situation we may observe in 
the theory of white dwarfs, see, for example, [10] and the 
references therein, as well as some additional comments 
in Section V, B. 

Unfortunately, we still do not have a complete set 
of observational data for a direct confrontation of the 
present-days relativistic theory with astrophysical obser- 
vations. In particular we do not have precise observa- 
tional data both for the mass and the radii of a given 
neutron star, we do not know the precise upper limit of 
neutron star masses, e.t.c. 

The existence of such universal upper limit is a basic 
prediction of the modern relativistic theory of stars, but 
the theory is not able to give a definite predictions for the 
corresponding value, due to the uncertainties in EOS. As 
a result, all observed massive compact dark objects with 
gravitational mass > (5 — 10) x m Q are automatically 
interpreted as a candidates for black holes, despite of the 
fact that there still do not exist undisputable direct ob- 
servational evidences for existence of such exotic objects 
with their non-avoidable attribute - the event horizon 
[11]. 

The fast development of this scientific domain calls for 
a further investigation of different aspects of the general 
theory, including its basic assumptions. 

In the present series of articles we will not consider 
the EOS problem, nor the complicated dynamical prob- 
lems of rotating stars, or stelar oscillations. Here we 
will reconsider some basic features of relativistic theory 
of gravity when applied to the study of stelar physics in 
the simplest static spherically symmetric case. We shall 
show that there exist new classes of models for relativis- 
tic stars, thus enlarging essentially the general theoretical 
scheme. We hope that the new relativistic models may 
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lead to a better understanding of the real astrophysical 
observations. 

Here we utilize a new approach to the spherically sym- 
metric static solutions to Einstein equations (EE), based 
on careful analysis of the radial gauge. Recently this ap- 
proach brought into the world a new two-parameter fam- 
ily of solutions of EE with a massive point source and 
some unexpected physical consequences, see [12]. One of 
them is that in the point particle problem the global ana- 
lytical properties of the solutions of EE in complex plain 
of the radial variable are fixing this variable in a unique 
way, together with the corresponding boundary condi- 
tions. Similar phenomenon is well known in the theory 
of analytical functions: they are unambiguously defined 
by their singular points in complex domain. 

In mathematical sense this way were derived the fun- 
damental static spherically symmetric solutions of EE. 
These solutions are analogous to the fundamental solu- 
tions of classical Poison equation with point source. Here 
we are extending this approach to the theory of relativis- 
tic stars. 



II. THE PHYSICAL CONSEQUENCES OF THE 
CHOICE OF RADIAL GAUGE IN THE STELAR 
PROBLEM 

The EE determine the solution of a given physical 
problem up to four arbitrary functions, i.e., up to a choice 
of coordinates. This reflects the well known fact that GR 
is a gauge theory. According to the standard textbooks 
[1], the fixing of the gauge in GR in a holonomic frame 
is represented by a proper choice of the quantities 

i=9^dx (V\9\9 Xv ) , (II.l) 



which emerge when one expresses the 4D d'Alembert op- 
erator in the form g^V„V v = (d p d u - f^A). Un- 
fortunately, up to now physically reasonable principles 
for the choice of the gauge in GR are not known. More- 
over, at present many of the relativists are thinking that 
this is not a physically essential GR problem. 

We shall call the change of the gauge fixing expressions 
(II.l), without any preliminary conditions on the analyt- 
ical behavior of the used functions, a gauge transforma- 
tions in a broad sense. This way we essentially expand the 
class of the gauge transformations, we intend to discuss, 
looking for a physically meaningful choice of the gauge 
conditions (II.l). 

In the static spherically symmetric problems the struc- 
ture of the space-time is M^ 3 ) = T t (l) x x SO(3). 
There exists unambiguous choice of the global time t on 
the ID time-translations group T t (l) and of the angle 
variables 6, (f> - on the 50(3) group space. These vari- 
ables are unambiguously fixed by symmetry reasons. In 
proper units (in which the velocity of light is c = 1) this 
choice yields the familiar form of the space-time interval: 



with unknown functions gtt{f) > 0, g r r{r) < 0, p(r). 

Thus the form of three of the gauge fixing coefficients 
(II.l): r t = 0, Fe=— cot#, = is fixed by symmetry 
reasons, but the quantity 
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(II.3) 



ds 2 = g tt {r)dt 2 + g rr (r) dr 2 - p(r) 2 (d9 2 + sin 2 9 d<j) 2 )(II.2) 



and, equivalently, the function p(r) are still not fixed. 
Here and further on, the prime denotes differentiation 
with respect to the variable r. 

The physical and the geometrical meaning of the radial 
coordinate r is not defined by symmetry reasons and is 
unknown a priori [12, 13]. The only clear thing is that its 
value r = corresponds to the center of 50(3) symmetry. 
In the case of relativistic stars with regular distribution 
of matter this 3D-space point is the physical center of the 
star, where the mass m(r), surrounded by a sphere with 
coordinate radius r = 0, is m(0) = 0. 

We shall use this physical property of the mass as a def- 
inition of the star's center C, because it does not depend 
on the choice of the radial variable r. Thus the mass m 
can be used to find the geometrical place at which the 
proper radial variable r must equals zero. This is the 
main specific feature of our approach to the theory of 
the relativistic stars. 

In contrast to the radial variable r, the quantity p has 
a clear geometrical and physical meaning: p defines the 
area A p — Airp 2 of a centered at the center C sphere 
with "area radius" p. From physical point of view one 
can refer to this quantity as "a luminosity variable " (or 
"a luminosity radius"), because the luminosity of distant 
physical objects is reciprocal to A p . In other words, the 
variable p describes the spherically symmetric spreading 
of energy of any kind. 

We refer to the choice of the function p(r) as a choice 
of radial gauge in a broad sense [12], allowing, in general, 
singular changes of the variable r. We call the freedom 
of choice of the function p(r) " a rho-gauge freedom" in a 
broad sense, and any definite choice of function p(r) - " a 
rho-gauge fixing" . 

At first glance the fixing of the function p(r) seems to 
be rather arbitrary and without any physical significance. 

From geometrical point of view the choice of the ra- 
dial gauge defines an imbedding of the ID quotient space 
m[ 1} = (M(^ 3 )/T t (l)) /50(3) into the space-time M^ 3 ). 

For fixing of this additional mathematical structure one 
needs some physical conditions like boundary conditions, 
or conditions for fixing of the number and the charac- 
ter of singular points of the solution of EE in the whole 
complex domain. This was demonstrated in [12] for the 
case of fundamental singular solutions of EE with mas- 
sive point source. These additional conditions play an es- 
sential role in the problem, because they are determining 
the global analytical properties of the solutions. Actually 
they define the very manifold M^ 1 ' 3 ). 

The EE are holomorphic ones and their solutions must 
be studied in the whole complex domain of correspond- 
ing variables. The very EE do determine only the local 
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structure of M^ 1,3 ) . In our case the change of the function 
p(r) will be not a simple change of the labels of space- 
time points, if it changes the additional conditions, which 
fix the analytical properties of the manifold M^ 1 ' 3 ) in the 
whole complex domain [12]. 

Our present consideration illustrates this important 
juncture on the more physical example of solar models: 

It is obvious that physical results of any theory must 
not depend on the choice of the variables. In particular, 
these results must be invariant under changes of coor- 
dinates. This requirement is a basic principle not only 
in GR. It is fulfilled for any already fixed mathematical 
problem. 

Nevertheless, the change of the interpretation of the 
variables may change the very mathematical problem and 
some physical results, because we are using the variables 
according to their interpretation. For example, if we are 
considering the luminosity variable p as a radial variable 
of the problem, it seems natural to put the center of the 
star at the point p = 0. In general, we may obtain a 
physically different stelar model, if we are considering 
another variable r as a radial one: in this case we shall 
place the star's center at a different geometrical point 
r = 0, which now seems to be the natural position of the 
physical center C . The relation between these two geo- 
metrical "points" and between the corresponding stelar 
models strongly depends on the choice of the function 
p(r), i.e. on the radial gauge. 

Thus, applying the same physical requirements, like 

Hat the center C = °' ( IL4 ) 

in different "natural" variables, we arrive at different 
physical models, because we are solving EE under dif- 
ferent boundary conditions. One has to find a theoret- 
ical or an experimental reasons to resolve this essential 
ambiguity, or one has to accept it as an non-avoidable 
component of the theory, recovering its physical meaning 
and its proper usage. 

The physical center C of the star is placed at the point 
r = by definition. To what value of the luminosity 
variable pc — p(0) corresponds the real position of the 
center C is not known a priori. This depends strongly 
on the choice of the rho-gauge function p(r). One can 
not exclude such nonstandard behavior of the physically 
reasonable gauge function p(r), which leads to some value 
PC > [12]. 

This very interesting novel possibility emerges in 
curved space-times due to their unusual geometrical 
properties and is not supported by our Euclidean ex- 
perience. It will be the main subject of study in the 
present series of articles. Such possibility was discovered 
at first in the original pioneering article by Schwarzschild 
[14] and discussed by Brillouin [15], but at present it is 
widely ignored. A physical necessity of considering val- 
ues of the p- variable, not less the Schwarzschild radius, 
was stressed by Dirac [16], too. 

The present-days standard theory of relativistic stars 
is based on the Hilbert radial gauge (HG): p(r) = r. In 



this gauge the center of the star is placed at the point 
pc = p(0) = 0. This rather arbitrary additional condi- 
tion was at first utilized by Schwarzschild in his simple 
model of incompressible stars [4] . There he had used for- 
mal mathematical reasons to be able to fix this way one 
of the integration constants. Actually he had postulated 
the global geometrical properties of the stelar center C in 
curved space-time, adopting the ones, which take place 
in the non-relativistic Euclidean case. 

The local reason seems to be Eq. (II. 4), which entails 
asymptotically flat 3D metric in a small enough vicinity 
of the stelar center C. According to GR, the spherically 
symmetric distribution of the masses outside this vicinity 
does not influence the flat geometry around C. 

Nevertheless, one has to take into account that there is 
no guaranty that starting from some luminosity ^ 
at the stelar surface, and going trough the curved 3D 
space back to the center C, defined by Eq. (II. 4), we 
will reach the value pc = 0. This is what we mean by 
" global" property of center C. This property of the cen- 
ter C depends on properties of the interior solution of 
stelar problem in the whole interior domain. From point 
of view of TOV system of differential equations, such 
property of the center C depends on the global proper- 
ties of the inner solution. 

The assumption pc — has a strong influence on the 
further development of the theory of relativistic stars. In 
particular, it forces one to impose a regularity condition 
at the point p = both on the matter distribution and 
on the solutions of EE. As a result, one uses only a very 
specific solutions of TOV equations for relativistic stars 
[l]-[8]. These solutions form a set of zero measure in the 
variety of all solutions of the problem. Thus one is forced 
to ignore the vast majority of the solutions, which are of 
general type and do not obey the regularity condition at 
the point p = 0. 

The novel solutions of EE for massive point particle, 
discovered in [12], raise a new understanding of the role 
of the luminosity variable p. Here we shall show that in 
the relativistic theory of stars the same approach allows 
a consideration of all solutions of the TOV equations in 
a physically meaningful way. This way we essentially 
enrich the relativistic theory of stars. 

An extremely important consequence of our more wide 
treatment of stelar models is the existence of relativistic 
stars with arbitrary large mass, and, at the same time, 
with arbitrary small geometric radius and arbitrary small 
luminosity. This unexpected possibility will be mathe- 
matically proved for incompressible relativistic stars in a 
subsequent article. 

The strong nonlinearity both of the differential equa- 
tions and the boundary conditions may yield, in gen- 
eral, several different classes of new solutions. This re- 
sembles the real situation, illustrated by the well known 
Hcrtzsprung- Russell diagram for stars in Nature [17]. 

Because of presence of the additional parameter pc in 
the general solutions of TOV equations, the total mass of 
the star to* is not a function only of its coordinate radius 



4 



r* (or geometrical radius i?*, or luminosity radius p*) and 
may vary independently of it. Due to this property, our 
approach gives for the first time a possibility to consider 
the point particles in GR as a limiting case of a body 
with finite dimension, much like in the Newton theory of 
gravity. This important new feature will be described in 
another subsequent article. 

The considerations in this first of series of articles has a 
preliminary character, setting in a new way the relativis- 
tic theory of gravity in the stelar physics. This article 
does not aim a construction of a specific models of rela- 
tivistic stars. Here we give only the new general scheme 
for such considerations. Due to technical reasons we de- 
scribe specific models of relativistic stars with different 
EOS and other further developments elsewhere. 



III. SOLUTION OF THE EXTENDED TOV 
SYSTEM IN HILBERT GAUGE 

The luminosity variable p gives a very convenient de- 
scription of stelar structure in real domain. The use of 
this variable ensures a local radial-gage-invariance of the 
approach to solution of this problem. Therefore, keeping 
the traditions, we shall work in this Section in HG. In our 
consideration the values p < pc have no physical mean- 
ing. The values p G [pc,P*] describe the inner domain 
and the values p G (p»,oo) correspond to the exterior 
vacuum domain outside the star. 

Then the inner metric (II. 2) for a static spherically 
symmetric star is defined by the metric components 



gu(p;p*,Pc) = e 2 ^r-e c \ 
l-2m(p; p. t ,pc)/p 



(III.l) 



Here and further on we are using units c = Gjy = 1- 

The mass m(p; p*,pc), surrounded by a sphere with a 
luminosity radius p, obeys the first TOV equation: 



dm 
dp 



Airep > 0, 



(III.2) 



supplemented by the boundary conditions 



m{pc;p*,Pc) = 0, m(p*;p*,p c ) = m(p*,Pc)-(III.3) 

The first condition is the definition of the physical center 
C of the star, placed at a position with unknown value of 
the luminosity variable pc > 0. The second one defines 
the total mass of the star m(p*, pc), obtained using the 
unknown value of the luminosity variable p* > pc at the 
edge of the star. 

As a consequence of Eq. (III.l) and (III. 2) we ob- 
tain for any admissible value of pc the relation 

-9pp(pc;p*,pc) = l- 

The pressure p(p; p*,pc) obeys the equation 



dp (p + e) (m + 4irp 3 p) 
dp p(p — 2m) 



where e(p;p*,p c ) is the energy density. The boundary 
conditions for this equation are: 

p(pc;p*,Pc) = Pc{p*,Pc), p(p*;p*,pc) = 0. (III.5) 

Here pc is the pressure at the star center C. The second 
condition defines the physical edge of the star. 

One has to extend the above TOV system adding the 
equation for the proper mass mo(p; p*, p c ) of the star in 
the sphere with luminosity radius p: 



dmo 
dp 



= 4Trep 2 y/-g pp > 0, 



(III.6) 



together with the boundary conditions 

m (p c ;p*,pc) = 0, m (p*;p*,p c ) = m (p*,p c ) .(III.7) 

and the equation for the gravitational potential 
<fi(p\ P*,Pc)- 



dip m + 4:irp 3 p 



>0, 



(III.8) 



<0, 



(III.4) 



dp p(p — 2m) 
together with the boundary conditions 

^(pcsp*, pc)^ ( pc(p*, pc), p(p*;p*, pc)=<p(p*, pc)-(ni.9) 

To have a closed system of mathematical equations 
one has to add the EOS. It can be defined in different 
equivalent forms. The most convenient for our general 
considerations is the following one: 

e = e(p). (111.10) 

It is useful to introduce, too, the quantity 

w=p/e. (III.ll) 

As a result of our consideration we see that the descrip- 
tion of stelar structure leads to a correct mathematical 
boundary problem with unknown ends pc and p*. 



IV. SOLUTION OF A CAUCHY PROBLEM AS 
A METHOD OF SOLUTION OF THE STELLAR 
BOUNDARY PROBLEM 

As we have seen in the previous Section, the stelar 
structure is determined by solution of the boundary prob- 
lem, described by Eq. (III.2)-(III.10). It is remarkable 
that in the simple case at hand the subsystem of differ- 
ential equations (III. 2), (III. 4) splits and can be solved 
independently of the other equations in ETOV system. 
Using the first of the conditions (III. 3) and (III. 5) as 
initial conditions for this subsystem, one obtains the so- 
lutions of the corresponding Cauchy problem in the form 

[1]: 

m = m(p;p c ,Pc), P = p(p; Pc , Pc) ■ (IV.l) 
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Using the already known function p(p; pc-,Pc)i one can 
solve the second of the Eq.(III.5), written in the form 
p(p*: PCiPc) = with respect to the quantity pc- Thus 
one obtains the pressure at the center C in the form 
Pc = Pc(p*,Pc)- Now, substituting this function in 
the known expressions m(p; pc,Pc) and p(p; pc,pc), and 
solving the corresponding integrals, one obtains the so- 
lution of the whole boundary problem, described in the 
previous Section, in the form: 

m(p;p*,p c ) = m(p;p c ,Pc(p*,Pc)), (TV. 2a) 

p(p;p*,pc) =p(p;pc,pc(p*,pc)), (TV. 2b) 

i \ a f P e (P\P*-,Pc)p 2 dp m . oA 

m (p; P*,Pc) = 4tt / - (TV. 2c) 



he V 1 - m (p;p*,pc)/p 
<p(p; p*,pc) = <pc(p*,pc) + 

<" p m(p; p*,pc) + 4np 3 p(p; p*,pc) 



J 



p(p - 2m(p; p., pc)) 



dp. (TV. 2d) 



Here 



e(p;p*,Pc) = H^*^))' *' e ^'*j ; (IV.3) 

[0, if p e (p*,oo) 

is obtained making use of Eq.(III.lO). 

The exterior solution in HG is well known: 



9tt(p; P*,Pc) : 



-1 2m(p.,p c )_ 

9pp(p;p*,pc) p 



It gives an interpretation of the constant m(p*, pc) as a 
Keplerian mass of the star. 

The inner solution (IV. 2) depends on the function 
<fic(p*i Pc) which can be determined using Birkhoff the- 
orem. The gravitational field at the edge of the star de- 
pends only on the total mass m(p*, pc) and matches the 
exterior vacuum solution, which is unique, up to choice of 
radial gauge. When written in HG, t he matching condi- 
tion gives ip(p*,p c ) = lng tt * = In \J\ - 2m(p*, pc)/p*- 
Then from the last equation (IV. 2d) one obtains 



<pc(p*,pc) = In \J\ - 2m(p^ 1 p c )l 'p* - 
rP * m(p;p* 7 p c ) + 4irp 3 p(p;p*,p c ) 



p(p-2m(p;p*,p c )) 



dp. (IV.5) 



As a result 



f(P- P*,Pc) = In yjl- 2m(p*,p c )/p* 
'* rn(p; p*,pc) + 4?rp 3 p(p; p*,p c ) 
p(p-2m(p;p*,p c )) 



dp (IV.6) 



and the inner solution depends only on the two parame- 
ters pc and p* with unknown values. 

The geometrical radius of the star i?„ = R(p*,pc) = 
Ipc V-9p P dp = Jo* V~9rrdr is 



R(p*,pc) 



pdp 



Pc y/p(p-2m(p; p»,pc)) 



(IV.7) 



Thus for a given EOS we arrived at a two-parameter 
family of relativistic stars. Up to now the solutions have 
been parameterized by luminosity variables pc and p*. 

A similar procedure, based on solution of back Cauchy 
problem with initial point at the edge of the star, p*, 
illustrates in the best way our physical definition of the 
center of star C: 

We can solve the subsystem of differential equations 
(III. 2), (III. 4) under much more physical initial condi- 
tions - fixing in arbitrary way the directly measurable 
mass > and luminosity radius p* > pc = 2m* 
and using the value of pressure p* = at the edge of 
the star. Now we can integrate the differential equa- 
tions back with respect to the variable p. According to 
Eq. (II. 4), the center C of the star is defined as a point 
pc < p*, at which m(pc, p*, pc) — 0. Finding this way 
Pc = Pc( m *, P*), we obtain the not-directly- measurable 
pressure p c (= Pc(p*,Pc) = p(pc, P*, Pc)) at the center 
C, which, itself, is hidden for us, from observational point 
of view. 

Obviously, the widespread in the literature [1] stiff re- 
lation m* = to* (p* ) will appear only if we pose by hands 
the commonly adopted extra condition pc(m*, p*) = 0, 
although there are no physical reasons to do this. 

It is clear, that the procedure, based on backward in- 
tegration, lies on much more physical ground, than the 
standard one. It is complete equivalent to the traditional 
procedure, if we use unknown value of the luminosity vari- 
able pc for solution of Cauchy problem with the center 
of the star C as starting point. 



V. THE MAPPINGS Rm AND Rm 

A. The Mapping Rm 

If one solves the algebraic equations 

R(p*,Pc) = R*, m(p*,pc) = m* (V.l) 

with respect to the variables p* and pc , expressing them 
as a functions of the variables i?* and m* , one obtains a 
complicated nonlinear mapping 



{i?*,m*} 



Rm 



{p*(R*,m*),pc(R*,m*)}. (V.2) 



The study of this mapping is a basic physical problem in 
our approach to the relativistic stelar structure. 

This way one can parameterize the solutions of ETOV 
equations for relativistic stars with a fixed EOS by the 
two parameters i?* and m*, which are directly measur- 
able. 

As we see, in our model of relativistic stars of most 
general type, the theory of gravity in HG does not yield 
a functional dependence between the mass to* and the 
radius i?*, only. For a fixed value of i?*, according to 
Eqs. (V.l), (V.2), the mass of the star to* can still vary. 

To obtain a stiff functional dependence between the 
mass to* and the radius i?* one must introduce some 
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auxiliary condition. In the commonly accepted models 
of relativistic stars the role of such condition plays the 
assumption pc = 0, which seems to be not necessary from 
physical point of view. If imposed, this extra condition, 
together with relations (V.2), gives the well known stiff 
functional relation to* = m*(.R*). 



B. The Mapping Rm 

If one solves the algebraic equations 

R(p*,Pc)=R*, m Q (p*, p c ) = to » (V.3) 

with respect to the variables p* and pc, expressing them 
as a functions of the variables i?* and too* , one obtains 
another complicated nonlinear mapping 

{i?* , too* } Rm °> {/?* (R* , to * ) , pc (R* , m-o* ) } • ( V.4) 

This way one can parameterize the solutions of ETOV 
equations for relativistic stars with a fixed EOS by the 
two basic parameters i?* and mo*. 

The study of the mapping Rm (V.4) is the second 
basic physical problem in our approach to the relativistic 
stelar structure in HG. 

In contrast to the Keplerian mass to*, which depends 
on the concentration of the fixed amount of matter in 
a given star, its proper mass too* is independent of this 
concentration and characterizes the very amount of mat- 
ter, the star is build of. 

If imposed, the extra condition pc = now yields a 
stiff functional dependence i?* =-R*(too*) for any given 
EOS. Combined with the relation m* =m*(i?*) from the 
previous Subsection, it leads to another stiff dependence: 
to* =to*(to *). 

The above stiff relations are the most important spe- 
cific prediction of the relativistic theory of gravity with 
auxiliary condition pc = 0. 

We have to stress that in Newton theory of gravity, 
which is known to describe well enough the physics of 
large class of real stars [1], including the Sun [18], as 
well as some features of white dwarfs [3], [10], we have 
analogous stiff functional relations, with similar origin - 
the regularity condition at the Euclidean point rc = 0. 

At the same time the functional dependance between 
corresponding quantities in the Newtonian theory is es- 
sentially different, in comparison with the standard rel- 
ativistic theory. In particular: a) There we do not have 
two different masses to* and mo, , because to* = too, ; 
b) In Newtonian models of spherically symmetric bodies 
with different EOS, as a role, we have no limitations on 
the mass to* and the radius i?*, due to the requirement 
to have a regular solutions at the stelar center C. The re- 
strictions on the mass to* appear, as an exception, only in 
the limit of degenerate ultra-relativistic matter [3] , which 
is not a realistic case. 



This is in a sharp contrast to the GR theory of stars, 
based on the condition pc = 0, in which we have restric- 
tions on the mass to* for any EOS [1], [3], [4], [19]. 

In our more general geometrical models of GR stars 
the regularity condition at center C with rc — are 
satisfied, because pc = p(0) ^ 0. Here we do not have 
stiff relations of the discussed type, without imposing 
some additional extra conditions. 

To check the existence of stiff relations between to*, 
too* and _R* in Nature, one has to analyze properly the 
observational data. The absence of stiff relations would 
lead to a dispersion of the observational data in a large 
domain of the corresponding variables. In the opposite 
case - if the stiff relations take place in physical reality, 
the observational data have to show a clear functional de- 
pendence between the corresponding quantities for some 
class of real stars with fixed EOS and matter content, 
which are in the same instant state. 

This phenomenon can help us to test the validity of 
relativistic theory of gravity with pc = (or with some 
other extra condition) in stelar physics, performing a pre- 
cise analysis of the observational data. 

Even a cursory look at the Hertzsprung-Russell dia- 
gram [17] will convince us that the observations may not 
support the standard relativistic theory of stars in HG 
with the extra condition pc = adopted: 

In the Hertzsprung-Russell diagram we see a big dis- 
persion of the temperature-luminosity positions of stars 
with different masses to* and radii i?*. Unfortunately it 
is not clear whether the (non)existence of stiff relations 
can be mask completely by the strong dependence on 
EOS and instant time state of the star, which changes es- 
sentially during the time evolution. On the Hertzsprung- 
Russell diagram we are witnessing some mixture of dif- 
ferent effects, due to too many physical factors. This is 
a serious obstacle for making some definite conclusions 
about the problem, we are discussing. 

Without any doubts, the best candidates for such anal- 
ysis are the white dwarfs. Their EOS is well known and 
fixed. In addition, some observational information for 
their radiuses and masses is available [22]. 

The corresponding stiff mass-radii relations were estab- 
lished on the basis of theory of degenerate stelar matter 
and studied in details in [3]. According to Provencal et 
al. (1998) [10], "One might assume that a theory as basic 
as stellar degeneracy rest on solid observational ground, 
yet this is not the case. Comparison between observation 
and theory has shown disturbing discrepancies 

Actually, a relatively large dispersion of observational 
data for masses and radii of white dwarfs are observed 
[10]. Its explanation, on the basis of standard relativistic 
theory of stars, forces one to accept a doubtful variations 
of the matter content of the white dwarfs. For example, 
a possible explanation of too small radii of some white 
dwarfs is the assumption about the existence of iron- 
reach core in them. According to Panei et al. (2000) 
[10]: "Obviously, such result is in strong contradiction 
with the standard predictions of stelar evolutionary cal- 
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culations, which allow for an iron-rich interior only in the 
case of prcsupernova objects" . 

Taking into account this situation, it seems interesting 
to analyze the existing data for white dwarfs mass-radii 
relation from point of view of the two parameter family of 
novel geometrical models of relativistic stars, presented 
here. We intend to perform such analysis elsewhere. 



The 2D Domains 



j(2) 

"P* ,PC 1 



4 2 > 



In our new model of relativistic stars the parame- 



ters {p*,pc} <E 
m>( 2 ) 



,(2) 

"P-,PC 



C 



vary in some 2D domain 



pj. pc , restricted by the conditions: 

< pc < P*, 
< 2m(p*,pc<) < p*, 

< Pc(p*,Pc) < oo. 



(V.5a) 
(V.5b) 
(V.5c) 



These conditions determine the physical 2D domains 



,(2) 

V R. ,m. 



and 



of the stelar parameters i?* > 0, 



to* > and to * > in the mappings (V.2) and (V.4). 



The form of the domains 



and 



])( 2 ) D (2) 

depends on the EOS. In general, its determination is a 
hard theoretical problem. Its solution is important for 
observational tests of the existence of the stiff relations, 
described in the previous two Subsections. 



VI. SCALE PROPERTIES OF THE ETOV AND 
HG SCALE INVARIANT QUANTITIES AND 
RELATIONS 

An important general property of the equations 
(III.2)-(III.8) was discovered by Bondi in 1964. This is 
their formal invariance under the scaling transformations 
with a constant coefficient A [6] : 

p^Ap, to^Ato, TO ^Am , £^A~ 2 e, p^\~ 2 p. (VI. 1) 

(See, too, the articles by Hartle, by Ellis et al., and by 
Collins in [6].) 

It is obvious that the quantities g tt , <p, g PP and w are 
A- invariant. 

If, and only if, w = const, the EOS (III. 10) is A- 
invariant and the solutions of the whole ETOV sys- 
tem will have a self-similar behavior, see the articles by 
Collins and by Rendal&Schmidt in [6]. 

Instead of the local binding energy Am(p, p*, pc) '■= 
too — to, which is not A -invariant, one can consider the 
ratio 

g{fi, />*, pc) :=m/m e (0,1). (VI.2) 

It measures in a A -invariant way the local mass defect 
of the star mater, i.e. the mass defect in the sphere with 
luminosity radius p and center C. 



Another important A-invariant local (in the above 
sense) quantity is f(p,p*,Pc) = g{p,P*,Pc) 2 - 
9tt{p, P*i Pc) + 1. In the case at hand it has the form 



f(p,P*,Pc) = ( — 
Vto 



+ 



2to 



= Q 2 +e, (vi.3) 



where ^ 2 = y > is the local compactness of the star. 

Considering the values of the corresponding quanti- 
ties at the edge of the star, one can introduce their 
global counterparts: £* := g(p*,pc) = g{p*,P*,pc), 
<;* := c(p*,Pc) = s(p*,P*,Pc) and /* := f(p*,pc) = 
f(p*,P*,Pc)- 

We are considering in details only the scale properties 
of the solutions of ETOV system for relativistic stars. 
The corresponding non-rclativistic equations have the 
same scale properties, because they can be considered 
as a special case of the relativistic ones, taking the limit 

c — > oo [1]. 



VII. THE SOLUTION OF THE ETOV SYSTEM 
IN BASIC REGULAR GAUGE 

A. The General Properties of BGR Inner Solution 

The basic regular gauge (BRG) is defined by the con- 
dition f r = 0. It has been proved to have a unique and 
important mathematical and physical properties [12]. 

Together with Eq. (H-3) the BRG definition gives 
a second order differential equation for the function 
Pbrg{t)i supplemented by the boundary conditions 



Pbkg(O) = pc, PBRc(r*) = p* 



(Vll.l) 



with unknown value of the radial variable r* . A simple 
integration of the differential equation gives: 



P 



-9pp 

gtt 



const. 



(VII.2) 



After one more integration of equation T r = 0, in the 
inner domain p <G [pc, P*] we obtain the relation: 



r 
r* 




(VII.3) 



We have used the boundary condition (VII. 1) at the cen- 
ter C and at the edge of the star to fix the unknown 
integration constants after the integration of Eq. (VII.2). 

The equation (VII.3) fixes the BRG function pbrg(t) 
in the interior of the star in the form 

pBR G (r) = Pbrg (jr'>P*>Pc) ' forre [0,r*],(VII.4) 

and yields the following basic properties of this function: 
i) Pbrg iv,P*,Pc) G [pc,P*] for n = e [0,1]; 
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ii) Pbrg (0;P*,Pc) = Pc\ 

iii) P 1 brg ( 1 ;p*.pc) = p*; 

iv ) PbAg (»7; PO) Po) = Po for rj e [0, 1]. 

v) dp/ dr7 > 0. 

vi) pc + {p*-pc)il<P l £ t RG (v,P*,Pc) <P*-(P»-Pc) (1-?7) S 




FIG. 1: The narrow constrained physical domain of the func- 
tion p = Psfl G (r);p,,pc), defined by its properties i) - vi). 

These general properties entail the representation: 
Pbrg {v,P*,Pc) = Pc+{p*-Pc)v+{p*-Pc)v(1-v)q(v,P*,Pc) 

with some nonnegative function q(r];p*,pc) € [0,1], 
which is bounded and continuous in the 3D domain 
{r?e [0,1], p. > PC >0}. 

Actually the function q(r];p*,pc) is to be find from 
Eq. (VII.3) and depends on the HG solution of ETOV 
system with given EOS. As seen in Fig.l, this dependence 
is quite weak. 

Replacing the luminosity variable p in corresponding 
HG-expressions (IV. 2) by the function (VII. 4), in the in- 
terior of the star we obtain the stelar quantities in BRG: 



m = m 



BRG 



IP*, PC , m = m 



BRG 



V = P 



BRG 



„BRG 



;p*,pc , £ = £ 



BRG 



r* 



-;p*,pc , 



-;p*,pc 



f'= I— ;p*,Pc) ; - for r G [0,r*]. (VII.5) 



B. The Matching of the BGR Inner and Exterior 
Solutions 



The BRG-solution in the exterior vacuum domain can 
be obtained making use of Birkhoff theorem in corre- 
sponding BRG radial variable r. According to [12] it is: 



PBR G (ry 



2m, 



l-Qie 2r/r 



for re [r„m.In(l/e.)).(Vn.6) 



Now we have to match interior and exterior solutions 
using proper physical requirements on the luminosity 
variable p: 

1) Because the luminosity of physical source of any 
kind of radiation is L = const / An p(r) 2 , we see that in 
absence of surface sources of the corresponding radiation 
the function p(r) must be continuous. Otherwise we will 
destroy the local energy conservation at the place, where 
p(r) has a jump. 

2) The jumps of the derivative p'(r) will induce jumps 
of the radial derivative of the luminosity, L'(r). This 
means a jumps of the radial derivative of surface energy 
density. Such phenomenon is physically possible only in 
presence of some surface agent, like surface force (surface 
pressure) due to some surface tension. Hence, in this case 
the star will have some thin crust. Similar phenomenon is 
familiar, for example, from theory of neutron stars [1, 8], 
where the crust is introduced and studied from different 
point of view. 

The presence of the crust will obviously yield an ob- 
servable consequences. Indeed, in this case, according to 
Eq. (VII. 2), the coefficient f r in d'Alembert operator 
9^>Nv = 9^ u (dpd v - Tpd,,) will be singular at the 
stelar surface: 



f r = -rs{r-n). 



(VII. 7) 



Here S(r) is ID Dirac delta function. We refer to the 
quantity F as to "crust parameter". Its role in field prop- 
agation through the stelar surface will be considered in 
the next Subsection. 

If we exclude the presence of stelar crust and corre- 
sponding jump of the radial derivative L'(r), the crust 
parameter will be zero: F — and the derivative p'{r) of 
the luminosity variable will be continuous function 

It is obvious that the justification of the matching con- 
ditions is impossible without right physical interpretation 
of the luminosity variable p. 

The above physical considerations entails the following 
mathematical consequences: 

1. The continuity condition PBRG(r*-0) =PBRG(r*+0) 
gives the basic relation 



r* = m* In 



1 



2771, 
P* 



(VII.8) 



The finite value := m* ln(l/p*)) corresponds to the 
physical space- infinity in BGR, i.e., to the geometric 
place of points in 3D space, where the luminosity vari- 
able p(r) — > oo for r — > r M — and the 4D spacetime is 
asymptotically flat. 

Hence, as in the case of point particle source of gravity 
[12], in BRG we have to consider the finite interval r G 
[O,?*^) as a real physical domain of the radial variable 
r. In comparison with the point particle problem, the 
difference is that in stelar models we have different forms 

(VII. 4) and (VII. 6) of the gauge function pbgr{ t ) in 
the interior domain of the star and in the exterior vacuum 
domain. 
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2. One can easily find that for the exterior solution 
(VII. 6) the value of the constant in Eq. (VII. 2) is const = 
I /ml. 

Taking into account that: 

a) The quantities p, g pp and g tt are continuous func- 
tions of the radial variable at the point r = r*; and 

b) The derivative p' BRG {r) > is positive everywhere 
in the physical domain of the BGR-radial variable r; 
we can describe the jump of this derivative by the formula 

p' B R G (r*+0) = e- F P ' BRG (n-0), F € (-co, oo). (VII.9) 

Obviously, for F ^ this formula describes a refraction 
of the lines p = p(r) at the point r*. It yields the relation 



mj 2 f p * dp 
r * J pc P 2 



(VII. 10) 



Now, making use of the already found HG functions 
m*( / o*,pc) an d m o(p*, Pc), an d matching conditions 
(VII. 8) and (VII. 10), we can solve the algebraic system 
of four equations 



/ to* 
\to * 



exp 



TO*(/9*,pc) - TO* = 0, 

m (p*,Pc) - to-o* = 0, 
2r* \ 2m* 



1 = 0, 



m 2 r p -dp LiPP_ ue F = Q (vim) 

J pc p V gu 

for six unknowns m,,mQ»,pc,/)»,r»,F with respect to 
the first four of them. Thus we arrive at a new form 
of our solutions for relativistic stelar models with given 
EOS: 



TO* 

Pi- 



rn 
: P 



BRG 



(r*,F), m * 



, BRG 



(r*,F), 



?* G (r*,F), pc = P^ G (r*,r).(VII.12) 



This representation sheds a new light on the physi- 
cal meaning of the two-parameter family of relativistic 
stars, obtained in present article: The constant parame- 
ter F defines the properties of stelar crust. For different 
values of this parameter we obtain relativistic stars with 
different crusts. 

After all, if we fix the value of the parameter F, we 
will obtain one parameter family of relativistic stars, pre- 
cisely as in Newton theory of stars and in the standard 
relativistic approach to stelar physics [1], but without 
extra condition pc = 0. 

For example, postulating continuity of the derivative 
Pbrg (t) a * t ne stelar edge r = r*, we obtain F = 0. 

The existence of one parameter family of relativistic 
stars with arbitrary fixed value of the parameter F be- 
comes possible just for the sake of matching conditions 
(VII.8) and (VII. 10). The condition (VII. 10) replaces the 
HG extra condition pc — and produces a new type of 
stiff relations in stelar physics. 

It is obvious that one can impose only one of these 
alternative extra conditions. A novel problem in stelar 
astrophysics is to verify which one of them, if any, takes 
place in Nature. 



C. Spreading of Waves and Static Fields Trough 
the Stelar Crust 

The physical agent, which brings into being the stelar 
crust, changes the space-time geometry in accord to con- 
dition (VII.9). Therefore the presence of the crust will 
influence the spreading of all possible physical wave fields: 
scalar, electromagnetic, gravitational, spinor, e.t.c. 

In this subsection we will present a preliminary investi- 
gation of the spreading of waves and static fields through 
the stclar crust. Our aim is to reach qualitative under- 
standing of possible role of the stelar crust for field's dy- 
namics and statics. Therefore we consider in proper ap- 
proximation only the simplest case of scalar spherically 
symmetric field $(£, r). To distinguish the effects, caused 
by the stelar crust, here we neglect the interaction of the 
wave fields with the stelar matter. The exact treatment 
of this issue is a complicated problem. Its consideration 
requires first to have a complete solution for some specific 
background stelar model. 

As a result of Eq. (VII. 7) one obtains for field $(i,r) 
in BRG the following wave equation: 

g^V^ = g u <Ptt+g rr <S>rr + F5(r - r* )* r = Q(VII. 13) 

Owing to the continuity of functions g tt (r) and g rr (r) at 
point r*, we can replace them in a small enough vicinity 
of the crust with their constant values <? (r*) and g rr {r*)- 
Then, changing the corresponding scales of time and ra- 
dial variables : t — > \/ 9tt{r*) t and r — > \J g rr {r*) r, and 
using the properties of Dirac 5-function, we arrive at the 
differential equation: 



<£>tt-<&rr + F S(r - r*)$ r * = 0, 



(VII. 14) 



where $ r * = $ r (t, r*) is the value of the first derivative 
<& r (t,r) at the stelar edge. 

The general solution of this equation can be repre- 
sented in the form of Fourier integral: 



$(i,r) = r 

J —I 



dwe iw *i?(r;r*,w). (VII.15) 



The amplitudes R(r;r*,u>) are described by the general 
solution 



R(r; r*, uj) — i?(r*,w) cos (w{r — r*)) + 

sin (uj{r — r*)) 



i? r (r* , w) (l + F AO (w(r - r* )) J ■ 



(VII. 16) 



of the second order ordinary differential equation: 

R rr + lu 2 R~ F5(r - r*)i? r (r*, u). (VII. 17) 

Here the arbitrary functions i?(r*,w) = i?(r*;r*,w) and 
i? r (r*, lu) = i? r (r*; r*, uj) appear as integration constants 
of Eq. (VII. 17) and for corresponding values of w present 
the values of the function R(r; r*, w) and its first deriva- 
tive R r (r; r*, oj) at the edge of the star; O(x) is the Heavi- 
side step function. For our purposes we have to regularize 
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this generalized function [20], i.e., we have to prescribe 
some definite value O(0) to this function at the point 
x = 0. Then A6(x) := 6(x) - 9(0). 
Now it becomes clear that: 

1. The physical role of the stelar crust is to pro- 
duce a jump in the sin-mode of the stationary waves 
(VII. 16). At the same time the cos-mode remains un- 
changed, crossing the crust. 

2. If F = 0, i.e., in absence of stelar crust, both modes 
spread trough the edge of the star as a completely free 
waves. 

3. Choosing R(r*,u>) = and proper special values of 
0(0), one can obtain solutions, which describe stationary 
waves only inside the star: 



R rns (r; r„, w) = -FRr(r*, w)e(-w(r-r„)) 



sin (u)(r— r*)) 



for 6(0) = 1 + 1/F, or only outside the star: 



R° uts {r;r*,u) = FR r {r*,u)@{u(r-r*)) 



sin [lo 



(r-r*)) 



- for 6(0) = 1/F. 

Thus we see that the parameter 6(0) plays the role of 
reflection coefficient for the sin-mode. For proper values 
of this coefficient we have a total inner, or total outer 
reflection of the sin-mode by the stelar crust. 

4. In the static limit oj — > one obtains from Eq. 
(VII.16): 

$(r) = $(r.)+$ r (r.)(l + rAe(r-r.))(r-r»).(VII.18) 

The last formula shows that crossing the stelar crust with 
F^0, the static field $(r) is a subject of refraction. 

The existence of the stelar crust with the above proper- 
ties is a new specific prediction of our models of relativis- 
tic stars. It may have important consequences not only 
for the stelar physics and needs further careful study. 



VIII. THE SOLUTION OF THE ETOV SYSTEM 
IN PHYSICAL REGULAR GAUGE 

The regular change of the rho-gauge, defined by the 
fractional linear mapping of the interval r £ [0, r^) onto 
the whole interval r £ [0, oo): 



r — > to* ln(l/£*) 



r + m*/ln(l/g*) 



(VIII.l) 



brings us to the physical regular gauge (PRG) [12]. There 
the radial variable r varies in the standard semi-infinite 
interval [0, r^). (Note that we are using the same nota- 
tions r, r*, ... for radial variables in different p-gauges.) 
Then in PRG we have 

M Spbrg(V,P*,Pc) , if r£ [0,r*]; 
PP * G(r) -\2m,(l-e 2 -)-\ ifr e[ ^,oo), (VIIL2) 



where now 

r(r* + m„/ln(l/g*)) 



V = 



r*(r + TO*/Ln(l/0*)) 



£ [0,1] for r £ [0,r*],(VIII.3) 



r, =to* 1/ln 




■1/ln 



> 0.(VIII.4) 



In Eq. (VIII. 2) we are using the modified Newton gravi- 
tational potential [12]: 



ip G (r; to*, to ») := - 



TO, 



r + m*/ln(^)- 



(VIII.5) 



From Eq. (VIII. 2) one easily obtains the important 
inequality for the stelar parameters, written in the fol- 
lowing two useful forms: 



2m, 2m, 

U-.= qI + — -<i e> p*>- — - 
p* i - ei 



(VIII.6) 



This is a more strong restriction on the domain *wp,, pc 
than the inequality (V.5b). Actually the inequality 
(VIII.6) is a direct consequence of matching condition 
p(r* — 0) = p(r* +0), written in the following BRG-form: 



2 - 



+ 



2m* 
p* 



2m* 



for r* > 0. 



IX. SOME CONCLUDING REMARKS 

It is clear that the new approach to stelar structure, 
developed in the present article, calls for revision of many 
of widely accepted features of the GR theory of stars. 
The changes are not based on the critics of the very GR, 
but on more deep understanding of its applications, and 
on solution of some open problems in this theory, like the 
physical justification of the choice of GR gauges. 

In particular, it is obvious that we must apply the 
Birkhoff theorem in PRG only in the interval p £ 
[p*,oo) ^ r £ [r*,oo), i.e. in the exterior vacuum do- 
main outside the star. As a result, in this domain all local 
GR effects like gravitational redshift, perihelion shift, de- 
flection of light rays, time-delay of signals, etc., are gauge 
invariant and will have their standard exact values. The 
PRG metric for this domain can be found in [12]. The 
differences between predictions of our general models of 
stars and the standard ones, based on the assumption 
pc = 0, can not be observed in the local gauge invariant 
gravitational phenomena, which take place in the outer 
vacuum domain, surrounding the stars. 

Hence, our most general geometrical models have an 
essential impact only on the theory of the interior of rel- 
ativistic stars, and on theory of spreading of different 
physical fields in stars, and around the stars. 

The obtained new results seems to deserve further 
study and can lead to serious changes of our understand- 
ing of physics of stars in Nature. Specific models of the 
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described novel type for relativistic stars with different 
EOS, as well as other developments and applications to 
problems of real stelar physics, will be published in sub- 
sequent articles of this series. 
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